Retarded interactions in Graphene systems 
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We first demonstrate how two-dimensional sheets are incorporated in the formalism for planar 
structures. Then we derive the interaction in the geometry of two free-standing graphene sheets 
and of one graphene sheet above a substrate. Numerical results are produced for the fully retarded 
interaction at K and at room temperature for undoped and doped graphene. Additional results 
are given both for a gold substrate and for an ideal metal substrate. 



I. INTRODUCTION 

The van der Waals (vdW) and Casimir interactions can 
be derived in many different ways. One way is to derive 
the interaction in terms of the electromagnetic normal 
modest of the system. For planar structures the interac- 
tion energy per unit area can be written as 
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dio 
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In [f k (iw)] , 



(1) 



where fk (uJk) = is the condition for electromagnetic 
normal modes. Eq. (|T|) is valid for zero temperature and 
the interaction energy is the internal energy. At finite 
temperature the interaction energy is Helmholtz' free en- 
ergy and can be written as 



S=I/ T |4E , ln[/ fe (ia; n )] 
w»=^5 n = 0, 1, 2, 



(2) 



where (3 = 1/fcgT. The integral over frequency has been 
replaced by a summation over discrete frequencies, so 
called Matsubara frequencies. The prime on the summa- 
tion sign indicates that the n = term should be divided 
by two. For planar structures the quantum number that 
characterizes the normal modes is k, the two-dimensional 
(2D) wave vector in the plane of the interfaces. The nor- 
mal modes can be found from self-sustained charge and 
current densities, from self-sustained potentials or from 
self-sustained fields. Here, we will find self-sustained 
fields. To do that one solves the Maxwell equations (ME) 
in all regions of the geometry and makes use of the stan- 
dard boundary conditions at all interfaces between the 
regions. For the present task we need a geometry con- 
sisting of three regions and two interfaces, 1|2|3. This 
geometry gives the mode condition function 



fk = 1 



e- 2 ^ M r 21 r 23 , 



(3) 



where ry is the amplitude reflection coefficient for a 
wave impinging on the interface between medium i and 
j from the i side, d is the thickness of region 2 and 



7i = yjl — Si (<j) (co/ck) . The function (oj) is the di- 
electric function of medium i and c is the speed of light 



in vacuum. For planar structures there are two types 
of mode, transverse magnetic (TM) or p-polarized and 
transverse electric (TE) or s-polarized. They have dif- 
ferent amplitude reflection coefficients. At an interface 
where there are no 2D sheet the TM and TE amplitude 
reflection coefficients are 



and 



TM 



TE 



(li-lj) 

hi + 7i) ' 



(4) 



(5) 



respectively. Note that r,-j = — ry holds for both mode 
types. 

In next section we show how these amplitude reflec- 
tion coefficients are modified at an interface containing a 
2D sheet. In Sec. Mil we discuss the dielectric function of 
undoped and doped graphene. Then follows the analyt- 
ical expressions and numerical results for free-standing 
graphene sheets, Sec. Sec. [Hi a graphene sheet above 
a gold substrate, Sec.|Vl and a graphene sheet above an 
ideal metal substrate, Sec. lVIl We end with summary and 
conclusions in Sec.lVIIl 



II. 



HOW TO INCORPORATE 2D SHEETS IN 
THE FORMALISM 



There are different formulations of electromagnetism 
in the literature. The difference lies in how the conduc- 
tion carriers are treated. In one formulation these car- 
riers are lumped together with the external charges to 
form the group of free charges. Then only the bound 
charges contribute to the screening. We want to be able 
to treat geometries with metallic regions. Then this for- 
mulation is not suitable. In the formulation we use the 
conduction carriers are treated on the same footing as 
the bound charges. Thus, both bound and conduction 
charges contribute to the dielectric function. In the two 
formalisms the E and B fields, the true fields, of course 
are the same. However, the auxiliary fields the D and H 
fields are different. To indicate that we use this alterna- 
tive formulation we put a tilde above the D and H fields 
and also above the dielectric functions. 
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With this formulation MEs have the form 

V • D = 4irp ext 
V-B = 

Vx ?+y| =° 4 



(6) 



and the boundary conditions at an interface between two 
media 1 and 2 are 



(E 2 - Ei) x n = 
(f) 2 - Di) ■ n = Mp s ) ext 
(B 2 -Bi)-n = 
(Ha-Hi) xn = -^K ei , 



(7) 



We note that the sources to the fields in MEs are the ex- 
ternal charge and current densities and also in the bound- 
ary conditions discontinuities in the normal component 
of the D fields and tangential component of the H fields 
are caused by external surface charge densities and ex- 
ternal surface current densities, respectively. The unit 
vector n is the surface normal pointing into region 2. 

The amplitude reflection coefficient gets modified if 
there is a 2D layer at the interface. We treat the 2D layer 
at the interface as external to our system. The fields will 
then induce external surface charge and current densi- 
ties at the interface. Two of the boundary conditions 
are enough to get the modified Fresnel coefficients. The 
other two are redundant. We choose the following two 

(E 2 - Ei) x n = 



(H 2 - Hi) 



xn = -% = ->x(Exn) 



(8) 

where a is the conductivity of the 2D sheet. The modified 
amplitude reflection coefficient for a TM mode is 



r TM _ £jH - £Hj + 2 7^7j" 11 



(9) 



where the polarizability of the 2D sheet is obtained from 
the dynamical conductivity 



a 11 (k,uj) = 



(k,co) k 



and the dielectric function is 



(k,u) = 1 + a 11 (k,uj) . 



(10) 



(11) 



For TM modes the tangential component of the electric 
field, which will induce the external current, is parallel to 
k, so the longitudinal 2D dielectric function of the layer 
enters. The bound charges in the 2D sheet also contribute 
to the dynamical conductivity and the polarizability. 

The modified amplitude reflection coefficient for a TE 
mode is 



r TE = 7»-7j +2(u>/ck) 2 a ± 
ij 7i+7j-2(w/cfc) 2 a- L ' 



(12) 



where the polarizability of the 2D sheet is obtained from 
the dynamical conductivity 



a 1 - (k,u) — 



2ma 1 - (k, lj) k 



and the dielectric function is 



e 1 - (k,u>) = 1 + a 1 - (k,uj) . 



(13) 



(14) 



For a TE wave the electric field is perpendicular to k, so 
the transverse 2D dielectric function of the layer enters. 
The bound charges in the 2D sheet also contribute to the 
dynamical conductivity and the polarizability. 

Now, when we have geometrical structures with 2D 
layers at some of the interfaces we just substitute the 
new reflection coefficients at the proper interfaces. Note 
that j*2 1 = — r\2 n0 longer applies so one has to be careful 
when making the substitutions. To arrive at the starting 
expression one might have used this relation. For the 
three layer system 1|2|3 we find the most general mode 
condition functions as 

f = i _ e -272faV 21 r 23 ; 
f™ = 



f TE = 



£l72-£27l+ 2 7l72<a^ 
£172 +£27i+ 2 7i72a^ 

72-7l+2(^/cfc) 2 q 



72+7l-2(w/cfe) 2 a^ 



£3 72~£ 2 73+27273«r 
£372 +£273 +27273 c*h 

72-73+2(^/cfc) 2 q^: 
72+73-2(w/cfc) 2 a^; 



(15) 

where we assumed there arc 2D sheets at both interfaces. 
These sheets may be different so we have put the sub- 
scripts L and R on the polarizability of the left and right 
sheet, respectively. 



III. DIELECTRIC FUNCTION OF GRAPHENE 

In terms of the polarizability the dielectric function 
of a 2D system is given by e (k, u) = 1 + a (k, to) 
1 



„2D 



(fc)x(k,w), where v 2D (k) = 2ire 2 /k is the 2D 
fourier transform of the coulomb potential and \ (k, w) 
the density-density correlation function or polarization 
bubble. It can be useful to have the dielectric function 
also expressed in terms of the dynamical conductivity, 
<r(k, u). The relation is: e (k, u>) = 1 + 2iri<r (k, u>) k/ui. 
These relations are valid for all 2D systems. 

From now on we limit the treatment to graphene. We 
discuss the T = OK functions, only, and argue later that 
we may use these functions also at 300K. All functions are 
longitudinal versions; transverse versions are not avail- 
able in the literature. All functions are derived under the 
assumption that the band structure for graphene consists 
of two pair of cones; each pair consists of a conical va- 
lence band with the tip in positive energy direction and a 
conical conduction band with the tip in the negative en- 
ergy direction; the cones are stacked upon each other so 
that the two tips coincide; the cones extend from minus 
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to plus infinity. In a real graphene system there are de- 
viations from the conical shapes and there are also other 
bands contributing. Thus the functions we discuss here 
are only strictly valid for small frequencies. We estimate 
the results to be reliable for Tiu < 4eV. 

Let us first begin with an undoped graphene sheet. In 
a general point, z, in the complex frequency plane, away 
from the real axis the density-density correlation function 



is 2 



X(k,z) = - 



16ft vW 



(16) 



where v is the carrier velocity which is a constant in 
graphene (E = ±hvk), and g represents the degeneracy 
parameter with the value of 4 (a factor of 2 for spin and 
a factor of 2 for the cone degeneracy.) In the numerical 
calculations we use the value— 8.73723 x 10 5 m/s for v. 

When the graphene sheet is doped the dielectric func- 
tion becomes much more complicated. However it has 
been derived by several groups^—. In the two next 
equations we use dimension-less variables: x — k/2kp; 
y = fou/2Ep; z = hz/2Ep. The density-density correla- 
tion function in a general point in the complex frequency 
plane, z, away from the real axis is— 



X (k, z) 



Ai{i 



4v ^F [*--/(M)]}; 

/ (x, ~z) = asin (! -£)+ asin( l±^) 
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(17) 
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where D = y gn/nfi 2 v 2 is the density of states at 
the fermi level and n is the doping concentration. The 
same result holds for excess of electrons and excess of 
holes. The expression in Eq. (|17[) reproduces the results 
of Refs. andi when z approaches the real axis from 
above, apart from an anomaly in the result of Ref.[H 
The expression in their Eq. (8), which should be purely 
real has an imaginary part. If this is removed their result 
agrees with both ours and that in Refill. 

In the calculations we need the function at the imagi- 
nary frequency axis. We have derived a very useful ana- 
lytical expression valid along the imaginary axis; an ex- 
pression in terms of real valued functions of real valued 
variables:— 



X' ( k ,w) = x(Mw) 
= -D \l 



4vV + 

g(x, y) = atan [h(x, y)k (x, y)] + I (x, y) 



h(x,y) = 



l{[x 2 { y 2 -l) + {y 2 + l) 2 ] 2 + {2yx 2 ) 2 } 



1/4 



\J {x 2 +y 


2 -i) 2 +(2y) 2 -(y 2 +i 
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sin |iatan 


2yx 2 




[x 2 ( y 2 -l) + (y 2 + l) 2 \ 


}■ 



y / -2x 2 (y 2 -l)-2(y4-6y 2 + l)+2(y 2 + l)y/~xi+2x 2 (y 2 -l) + (y 2 + l) 2 



(18) 



where the arcus tangens function is taken from the 
branch where < atan < n. The density-density cor- 
relation function on the imaginary frequency axis has 
been derived before in a compact and inexplicit form 
(see Ref.|6| and references therein.) Here we have chosen 
to express it in an explicit form in terms of real valued 
functions of real valued variables. 



IV. FREE-STANDING GRAPHENE SHEETS 

Two free-standing graphene sheets we obtain from our 
geometry by letting all three media be vacuum and let 
both sheets have the same polarizability. We further as- 
sume that there is negligible difference between trans- 
verse and longitudinal screening in graphene. Then the 
mode condition functions in Eq. (1151) reduce into 



jTM _ ^ _ e -27 <0) fcd 
jTE _ i _ e -2y ia) kd 



7 (0) q(fc,Lj) 



1-HyW a(k,u) 
(uj I ck) 2 a(k ,u) 



7(°>-(w/cfc) 



where 
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(0) 



(fc,w) = \Jl-(u/ck) 2 . 



(19) 



(20) 



If retardation effects are neglected there are no TE 
modes and the mode condition function for TM modes is 
reduced further into 



/ 



TM 



1 



-2kd 



a (k, uj) 
1 + a (k, uj) 



(21) 



Now we have all we need for writing down the expres- 
sion for the interaction energy between two free-standing 
graphene sheets. It is 



E 



-2-r m kd 



-f i0 >a(k,iuj) 



1+7(0) a(fc,iu) 



00 r 
^#/l^|l- e - 27(0,ferf 



— (uj/ck) a(k,iuj) 
j( )-\-(uj/ck) 2 a(k.iLo) 



in the full retarded treatment. To save space we have 
suppressed the arguments of the 7^ (k,iw) function. In 
the non-retarded treatment the energy reduces into 



E 



d 2 k 



dco 
2^ 



ln<M 



-2kd 



a (fc, iu) 
1 + a (k, iu) 



. (23) 

At finite temperature the frequency integration in 
Eqs.([22]) and (|23|) is replaced by a discrete frequency 
summation as described in Eq.([5]). 



A. Undoped graphene sheets 

Let us first discuss two undoped graphene sheets. With 
the particular screening in graphene it turns out that 
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FIG. 1. (Color online) The attractive non-retarded^ interac- 
tion energy between two graphene sheets. The lower straight 
line is for undoped sheets, while the bent curves are for doping 
densities 1 x 10 10 , 1 x 10 u ,l x 10 12 ,andl x 10 13 cm" 2 , respec- 
tively, counted from below. The upper thick straight line is 
the Casimir classical result for two ideal metal half-spaces. 




FIG. 2. (Color online) The attractive retarded interaction 
energy between two undoped graphene sheets. The dashed 
straight line is the TE contribution while the filled circles are 
the TM contribution. The lower solid straight line is the total 
contribution, clearly dominated by the TM contribution. The 
upper thick straight line is the Casimir classical result for two 
ideal metal half-spaces. 



a (k/X, ilu/X) — a (fc, iui) and the separation dependence 
of the non-retarded interaction in Eq. (|23[) becomes very 
simple. A change in dummy variables removes the only 
d in the integrand and produces the factor of d~ 3 in 
front of the integral. The result is shown as the lower 
thin straight line in Fig. [I] The upper thick straight 
line is the Casimir classical result for two ideal metal 
half-spaces. The remaining curves in Fig. [I] are the re- 
sults from doped graphene sheets with doping concen- 
trations 10 10 , 10 11 , 10 12 , and 10 13 cm -2 , respectively, 
counting from below. These curves approach asymptotes 
with ~ d~ 5 / 2 for large separations. These asymptotes 
will cross the ideal metal half-space result and in the 
doping case there are bound to be retardation effects. 
It turns out that also 7^ (k/X, iw/X) — 7^ (k, w) so 
the separation dependence of the retarded interaction in 
Eq. (|22p becomes equally simple. The same dummy vari- 
able change removes the only d in both integrands and 
produces the factor of d~ 3 in front of the integrals. In 
the undoped case the retardation effects are negligible. 
There is a small reduction in strength of less than one 
percent. The retarded results are shown in Fig. [51 Note 
that here the figure begins at 100 A. The TE contribution 
is negligible. If plotted in this hgure the non-retarded re- 
sult would fall within the thickness of the straight line 
representing the total result. Note that all four curves 
follow the same power law. 

At finite temperature we can no longer use any dummy 
variable substitution in the frequency variable, so we will 
not have a simple power law for all distances. Further- 
more we need the polarizabilities at finite temperature. 
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FIG. 3. (Color online) Same as Fig.fjbut now at 300 K. 



The change in polarizability with temperature is largest 
for small frequencies. At 300K the static value increases 
with approximately 10 percent 7 . The TE amplitude re- 
flection coefficient vanishes at zero frequency and the 
TM amplitude reflection coefficient increases by approx- 
imately one percent due to the finite temperature. The 
finite temperature correction is dominated by the n = 
term in the Matsubara summation that only depends 
on the static value. Thus we may safely use the zero 
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FIG. 4. (Color online) Same as Fig.[2] but now for doped FIG. 5. (Color online) Expanded view of Fig.g] where also 
sheets. the non-retarded result has been added (the thin solid curve 

with open circles.) 



temperature form of the polarizability. The results are 
shown in Fig. [3] The TE contribution becomes even less 
important at 300K. The TM contribution approaches a 
new asymptote ~ dr 2 . Also the Casimir classical ideal 
metal half-space result follows an asymptote with this 
slope at large enough separations. Note that the asymp- 
tote for the total result is somewhat smaller than half 
the ideal metal asymptote. It would be one half of the 
ideal metal result if the static polarizability were infinite. 
Using the approximation of a constant conductivity as in 
Ref.@ would give that result. 



B. Doped graphene sheets 

When the graphene sheets are doped the dielectric 
function becomes much more complicated (see Eq. (|18p .) 
It is no longer possible to use the simple change of dummy 
variables that we did before to find a simple power law for 
the interaction. We limit the presentation to the highest 
of the four doping concentrations in Fig. [I] The results 
at zero temperature are shown in Fig.[4j Here we see 
that neither the TM nor the TE contribution are straight 
lines parallel with the ideal metal result. Both increases 
so that the total result approaches the ideal case. To be 
noted is that the TE contribution now is as important 
as the TM contribution. They are responsible for half 
the result each at large separations. The total result 
does not cross the ideal metal half-space result; it stays 
below but close. This is due to retardation effects. Here 
it is interesting to compare to the non-retarded results. 
We show this comparison in an expanded part of the fig- 
ure as Fig.[S] Here we see even clearer that the TE and 
TM contribute equal amounts. The thin solid curve with 
open circles is the non-retarded result. It follows the to- 
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T = 300 K 
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FIG. 6. (Color online) Same as Fig.g]but now at T = 300 K. 



tal retarded result at the left end of the figure. Then it 
detaches itself and crosses the classical Casimir result. In 
this region the retardation effects are clearly present. 

At T = 300 K the TE result decreases in strength and 
vanishes for large separations. The results are shown in 
Fig. [6] The TM contribution then dominates completely 
for large separations. The TM waves result in a contribu- 
tion half the strength of the total contribution in the ideal 
Casimir geometry. This is because the static polarizabil- 
ity now diverges resulting in a TM reflection coefficient 
equal to unity. This behavior is completely analogous to 
the result in the real metal half-spaces geometry when 
the experimentally obtained dielectric functions are used 
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or the Drude model. Then the absence of TE contri- 
bution is due to dissipation. Here it is due to spatial 
dispersion. One of the benefits of geometries consisting 
of 2D sheets from the fundamental physics point of view 
is that spatial dispersion (momentum dependence of the 
dielectric function) is automatically fully included. To 
include spatial dispersion in a system of thick layers is 
very complicated—. In the large separation limit only the 
classical (n = 0) result survives. According to the Bohr 
- van Leeuwen theorem - there should be no classical TE 
contribution. Thus the theorem is fulfilled also in the 
present system. 



V. GRAPHENE SHEET ABOVE A GOLD 
SUBSTRATE 

A free-standing graphene sheet above a substrate we 
obtain from our geometry by letting media 1 and 2 be 
vacuum and medium 3 be the substrate. We further as- 
sume that there is negligible difference between trans- 
verse and longitudinal screening in graphene. Then the 
mode condition functions in Eq. (1151) reduce into 



fTM _ ^ _ e -27 (0) fcd 
jTE _ ^ _ e -2j m kd 



(u j ck) 2 a(k,u) 



s (w)7< >+ 7s 



7(°)-(w/c/c) 2 a(fc,ij) 



(24) 

If retardation effects arc neglected there are no TE modes 
and the mode condition function for TM modes is re- 
duced further into 



TM 



1 



-Ikd 



a (k, lo) 



1 + a (k, uj) 



e s (w) - 1 



e s (w) + 1 



(25) 



Now we have all we need for writing down the expression 
for the interaction energy for a free-standing graphene 
sheet above a substrate. It is 



E = h f yf4 f ^ In (l 



e -2 7 (0 »fcd 



~f (0 >a(k,im) 



,.00 



a(k,iu>) 



_ p -2 7 <°»fcd 



— (uj / ck) 2 a(k ,iu ) 
■y(°1+(uj/ck) 2 a(k,iuj) 
,(0 



X 



11 



<y(0) 



(26) 

in the full retarded treatment. To save space we have 
suppressed the arguments of the 7^ (fc, iuj) and j s (k, iw) 
functions. In the non-retarded treatment the energy re- 
duces into 



E 



2 Ay/ 



l-\-a.{k.iuj) 



s s (iw)— 1 
e s (iw)+l 



(27) 



At finite temperature the frequency integration in 
Eqs.([2U) and (|2"7|> is replaced by a discrete frequency 
summation as described in Eq.([5]). 
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FIG. 7. (Color online) The polarizability for gold at the imag- 
inary axis used in the numerical calculations. 
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FIG. 8. (Color online) The attractive non-retarded^ interac- 
tion energy between a graphene sheet and a gold substrate. 
The lowest curve is for an undoped graphene sheet while the 
other curves are for doping densities 1 x 10 10 , 1 x 10 11 , 1 x 
10 12 ,andlxl0 13 cm , respectively, counted from below. The 
upper thick straight line is the Casimir classical result for two 
ideal metal half-spaces. 



Before we continue we need to discuss the gold dielec- 
tric function, e s (icj). We have used experimental optical 
data— for the imaginary part of the dielectric function. 
We have extrapolated the data both for large and small 
frequencies and used a version 1 of the Kramers Kronig 
dispersion relations to find the result at the imaginary 
axis. The resulting data are shown in Fig. [7] 
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Undoped Graphene sheet 
above a gold substrate 




Undoped Graphene sheet 
above a sold substrate 



FIG. 9. (Color online) The attractive retarded interaction 
energy between a graphene sheet and a gold substrate. The 
dashed curve is the TE contribution and the filled circles are 
the TM contribution. The lower solid curve is the total contri- 
bution. The upper thick straight line is the Casimir classical 
result for two ideal metal half-spaces. 



A. Undoped graphene sheet 

The non-retarded results are shown in Fig.[8l The low- 
est thin curve is the result for an undoped graphene sheet 
above a gold substrate. We note that the result no longer 
follows the simple power law that we found for two free- 
standing graphene sheets. However, for large separations 
the curve approaches an asymptote ~ c?~ 3 and is par- 
allel to the classical Casimir result (upper thick straight 
line) . Included in the figure are also the results for doped 
graphene sheets. These curves are for doping densities 
1 x 10 10 , 1 x 10 11 , 1 x 10 12 , andl x 10 13 cm~ 2 , respectively, 
counted from below. The curve for the undoped graphene 
sheet will never come close to the classical Casimir result 
and we do not expect to find any important retardation 
effects. However, the results for the doped sheets will 
sooner or later cross the classical result and there should 
be important retardation effects. 

In the retarded treatment we see from Eq. (I2fj|) that in 
general it is not possible to make the dummy variable 
change that removes the only d in the integrand and pro- 
duces the factor of d~ 3 in front of the integral. However 
for very large d values the exponential factor in the log- 
arithm forces the frequency (and also k) to be small in 
order for the integrand to be non- negligible. This means 
that for large d values we may use the static value of 
the substrate dielectric function. When we do that we 
may use the dummy variable change and both the TM 
and TE contributions follow the power law d~ 3 . Thus 
for general separations the result will not follow a simple 
power law, but for large separations the result approaches 
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FIG. 10. (Color online) Same as Fig.Hbut now for 300 K 
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FIG. 11. (Color online) Same as Fig. [9] but now for a doped 
sheet 



an asymptote with the power law d~ 3 . This is what we 
find in Fig.O where we show the retarded results for an 
undoped graphene sheet above a gold substrate. Just as 
in the case of two free-standing undoped graphene sheets 
the TE contribution is much weaker than the TM contri- 
bution. On this plot the non-retarded results would fall 
on the curve for the total contribution. 

At 300 K (See Fig-HQ) the TE contribution becomes 
even less important. The TM contribution approaches a 
new asymptote ~ dr 2 . Also the Casimir classical ideal 
metal half-space result follows an asymptote with this 
slope at large enough separations. 
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FIG. 12. (Color online) Expanded view of Fig. [TT] where also FIG. 13. (Color online) Same as Fig.QT] but now for T = 
the non-retarded result has been added (the thin solid curve 300 K 
with open circles.) 



B. Doped graphene sheet 

Just as in the case of two free-standing doped graphene 
sheets the TE contribution becomes important and pro- 
duces half of the result for large separations. The results 
are shown in Fig.QTJ The total result does not cross the 
ideal metal half-space result; it stays below but close. 
This is due to retardation effects. Here it is interesting 
to compare to the non-retarded results. We show this 
comparison in an expanded part of the figure in Fig. 1121 
Here we see even clearer that the TE and TM modes con- 
tribute equal amounts. The thin solid curve with open 
circles is the non- retarded result. It follows the total 
retarded result at the left end of the figure. Then it de- 
taches itself and crosses the classical Casimir result. In 
this region the retardation effects are clearly present. 

At T = 300 K the TE result decreases in strength and 
vanishes for large separations. The results are shown in 
Fig. [13] The TM contribution then dominates completely 
for large separations. The TM waves result in a contri- 
bution with half the strength of the total contribution 
in the ideal Casimir geometry. The reason is that the 
static values of the graphene polarizability and the gold 
dielectric function diverges; the TM reflection coefficients 
for a vacuum-graphene interface and a vacuum-gold in- 
terface are both equal to unity; both TE reflection coef- 
ficients are zero because uj 2 a (k, iui) and uj 2 s s (iu)) go to 
zero when the frequency goes towards zero. 

This is completely analogous both to the result for two 
doped free-standing graphene sheets and to the result in 
the real metal half-spaces geometry when the experimen- 
tally obtained dielectric functions are used or the Drude 
model. Here it is due both to spatial dispersion and to 
dissipation. Also in the present geometry the Bohr - van 



Leeuwen theorem is obviously fulfilled. 



VI. GRAPHENE SHEET ABOVE AN IDEAL 
METAL SUBSTRATE 

An ideal metal substrate can be used as a reference 
system. It gives an upper limit for the interaction with 
real substrates. It is much easier to use than a real sub- 
strate since one avoids the complicated derivation of the 
dielectric function of the substrate on the imaginary fre- 
quency axis using Kramers Kronig dispersion relations. 
We include an ideal metal substrate for these reasons and 
for an additional reason. We found^-, v.i.z., that the non- 
retarded result at least for an undoped graphene sheet 
above an ideal metal substrate does not exist. The inte- 
grals do no not converge. We here show that there are 
no convergence problems in the retarded treatment. 

The mode condition function for a free-standing 
graphene sheet above an ideal metal substrate we obtain 
from the results for a real substrate in Eq. (I%4"|) by let- 
ting the dielectric function of the substrate go to infinity. 
Then the TM and TE amplitude reflections coefficients 
for the vacuum-substrate interface become unity and mi- 
nus unity, respectively. We find 



f™ 

f TE z 



1 



1 



-y i0 >a(k,ui) 



— (uj/ck) 2 a(k.uj) 



7< )-(w/cfc) 1! a(fc,w) 



(28) 



If retardation effects arc neglected there are no TE modes 
and the mode condition function for TM modes is re- 
duced further into 



f™ = l-e 



-2kd 



a (k, lu) 



l + a(k,iv) 



(29) 
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FIG. 14. (Color online) Same as Fig. [9] but now for an ideal 
metal substrate 



FIG. 15. (Color online) Same as Fig. 1101 but now for an ideal 
metal substrate 



Now we have all we need for writing down the ex- 
pression for the interaction energy for a free-standing 
graphene sheet above an ideal metal substrate. It is 



E = h f f ^ m (l - e-^ ( ° )kd 

k 



7 (0) q(fc,^) 



e -2 7 <°»fcd 



l+7< )ct(fc,ja;) 
(oj / ck) 2 a(k 



■y<-°> + (u]/ckya(k,iiu) 



in the full retarded treatment. To save space we have 
suppressed the arguments of the 7 < - ^ (Jk,iw) function. In 
the non-retarded treatment the energy reduces into 



E 



d 2 k 
(2nf 



^InO 

2tt 



-2kd 



a (k, ioS) 
1 + a (k, ioS) 



(31) 

This integral does not converge, at least not in the un- 
doped case. The reason is that the polarizabilty decays 
too slowly with frequency. 

At finite temperature the frequency integration in 
Eqs.([30|) and ([3l"j) is replaced by a discrete frequency 
summation as described in Eq.([5]). 



A. Undoped graphene sheet 

In the case of an undoped graphene sheet above an 
ideal metal substrate we see in Eq. ()30j) that we may use 
the dummy variable change in both integrals and both 
the TM and TE contributions follow the power law dr z . 
This is verified in Fig.[14]where the result for an undoped 
graphene sheet above an ideal metal substrate is shown. 
As can be seen the results are very similar to those for an 
undoped graphene sheet above a gold substrate, Fig.|Hl 
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FIG. 16. (Color online) Same as Fg.[TT]but now for an ideal 
metal substrate 



Here we can not compare to the non-retarded result since 
it does not exist. 

The results for 300 K are shown in Fig.Q~51 Again the 
results are very close to those for a graphene sheet above a 
gold substrate, Fig.JTU] However, for smaller separations 
outside the figures the total result for the gold substrate 
has a smaller slope. 



B. Doped graphene sheet 

The result for a doped graphene sheet above an ideal 
metal substrate is shown in Fig.[T5] As can be noted the 
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FIG. 17. (Color online) Same as Fg. 1131 but now for an ideal 
metal substrate 



results are very close to those for a doped graphene sheet 
above a gold substrate, Fig.HTl 

The room temperature results, T = 300 K, are shown in 
Fig.fTfl Again the results are very similar to the ones for 
a doped graphene sheet above a gold substrate, Fig. [13] 
However, for smaller separations outside the figures also 
here the total result for the gold substrate has a smaller 
slope. 



VII. SUMMARY AND CONCLUSIONS 

We have here demonstrated how to incorporate 2D 
sheets in the formalism for planar structures. We ap- 
plied the results to graphene systems. We found that 
for two free-standing undoped graphene sheets the retar- 
dation effects are negligible. The same holds for an un- 
doped graphene sheet above a gold substrate. However, 
for doped graphene sheets, two free-standing or one free- 
standing above a gold substrate, the retardation effects 
are bound to be important at large enough separations. 
However, at T = 300 K the retardation effects again go 
away. The retarded results for a graphene sheet above 
an ideal metal substrate are very similar to those for a 
graphene sheet above a gold substrate. On the other 
hand, the non-retarded results do not exist, so there we 
have a drastic difference. We should furthermore note 
that the "saturation effect" at small separations found 
in the interaction between two 2D metallic sheets^ is 
absent in geometries involving undoped graphene sheets. 
This spatial dispersion effect, present also in the case 
of doped graphene sheets in the region where the dop- 
ing carriers contribute, is due to the weakening of the 
screening in the sheet for large momenta. This happens 
for a momentum where the dispersion curve of the col- 



lective mode enters the single particle continuum. In the 
calculations for T = 300 K we used the zero tempera- 
ture polarizabilities for graphene. We estimate this to 
underestimate the large d asymptote with less than one 
percent. Furthermore, we used the longitudinal instead 
of transverse expressions for the graphene polarizabili- 
ties in the TE contributions; the reason for this is that 
nobody has derived the transverse versions, as far as we 
know; we plan to do this in the future. However, the TE 
contributions were found to be negligible in most situa- 
tions. 

Retarded results for graphene systems have been ob- 
tained by several group a 7 ' 12 " — before and Santos^ has 
discussed the unimportance of retardation in graphene 
systems. We will here go through these results one by 
one and compare them to ours. 

The interaction between an undoped graphene sheet 
and an ideal metal substrate was studied in Ref.[l2|. In 
their model they let graphene have a small but finite 
band gap and the carriers have non-zero mass. Thus, the 
comparison with our results is meaningless. 

In Ref.0 one studied the force between undoped 
graphene sheets at zero temperature. They used a model 
where the conductivity is a constant, <r(k, of) — <7q — 
e 2 /Ah. This is the limiting form for vk << w, 



a (k, iuj) = uja (k, itS) /2Trk = —uje 2 x (k, iuS) / k 2 
e 2 /4h, vk << ck < u>. 



me g k 



(32) 



One has to be careful with this approximation, since 
it is not valid in the finite temperature n = term, 
which dominates at large separations. They gave in their 
Eq. (36), for the first term of a perturbation expansion 
whereupon <7o is small, the expression E = —e 2 / (327rd 3 ) 
which is very close to our result. Note that it does 
not contain any material parameters. Furthermore they 
showed in their Fig. 6 how their conductivity varied with 
temperature; based on this result we concluded that we, 
to a good approximation, could use the zero temperature 
expression for the graphene polarizability in our room 
temperature calculations. 

In Ref.[l3| one studied the interaction between an un- 
doped graphene sheet and a substrate (gold and ideal 
metal) at finite temperature. They used a different ap- 
proach to the problem starting from a quantum-field- 
theory model at finite temperature and used different 
expansions. Our results show an overall agreement. How- 
ever, in the high temperature limit they put v = 0, which 
means that they use in that limit the same constant con- 
ductivity approximation used in Ref.0- This means that 
the graphene polarizability diverges at zero frequency 
and the TM reflection coefficient becomes unity. The 
result from this is that the interaction in this limit is 
half the corresponding value between two ideal metal half 
spaces, just as in the case of a doped graphene sheet 
above a gold or ideal metal substrate. We find a some- 
what smaller value (see Figs.fTUl and IT51) 

In Ref.[3 one studied the interaction between undoped 
and doped graphene sheets at zero temperature and at 
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300 K; also the interaction between a graphene sheet and 
a silica substrate was determined. Our results disagree 
with theirs. For two undoped graphene sheets they find 
many distance regions, each with a different power law. 
In the doped case the effect of doping sets in at very 
small separations with a fractional power law while in our 
results it sets in at much larger separations; their high 
temperature asymptote is much weaker than the ideal 
metal asymptote. We are convinced that their results 
are wrong. From the first appearance it looks as if they 
use exactly the same dielectric functions that we do for 
the undoped and doped graphene. After reading more 
carefully we believe that they have treated the strictly 2D 
sheets as layers of finite thickness, using the 2D dielectric 
functions in place of the bulk dielectric functions of the 



materials making up the layers. Then they have let the 
layers have a small but finite thickness. This procedure 
is not allowed. 

Finally, in Ref.[l5] the focus was on the change in force 
between a silica membrane and a metal plate when the 
membrane was covered by a graphene sheet. This work 
is just vaguely related to ours. We did not find anything 
that contradicts the findings in our present work. 
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